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Abstract. Dynamical chiral symmetry breaking (DCSB) is a remarkably effective mass generating 
mechanism. It is also, amongst other things, the foundation for a successful application of chiral 
effective field theories, the origin of constituent-quark masses, and intimately connected with 
confinement in QCD. Using the Dyson-Schwinger equations (DSEs), we explain the origin and 
nature of DCSB, and elucidate some of its consequences, e.g.: a model-independent result for 
OO [ the pion susceptibility; the generation of a quark anomalous chromomagnetic moment, which 

may explain the longstanding puzzle of the ai-p mass splitting; its impact on the behaviour of 
the electromagnetic pion form factor - thereby illustrating how data can be used to chart the 
momentum-dependence of the dressed-quark mass function; in the form of the pion and kaon 
valence-quark parton distribution functions, and the relation between them; and aspects of the 
O ' neutron's electromagnetic form factors, in particular F" /Ff and G"^. We argue that in solving QCD, 

^ I a constructive feedback between theory and extant and forthcoming experiments will most rapidly 

C ■ enable constraints to be placed on the infrared behaviour of QCD's j3 -function, the nonperturbative 

quantity at the core of hadron physics; and emphasise throughout the role played by confrontation 
, with data as a means of verifying our understanding of Nature. 

Keywords: Bethe-Salpeter equations, confinement, dynamical chiral symmetry breaking, Dyson- 
Schwinger equations, hadron form factors; hadron spectrum; parton distribution functions 
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1. INTRODUCTION 



A hundred years and more of fundamental research in atomic and nuclear physics has 
shown us that all matter is corpuscular, with the atoms that comprise us, themselves 
^ . containing a dense nuclear core. This core is composed of protons and neutrons, referred 

^ I to collectively as nucleons, which are members of a broader class of femtometre- scale 

particles, called hadrons. In working toward an understanding of hadrons, we have 
discovered that they are complicated bound-states of quarks and gluons, which are 
elementary and pointlike excitations whose interactions are described by a Poincare 
invariant quantum non-Abelian gauge field theory; namely, quantum chromodynamics 
(QCD). The goal of hadron physics is the provision of a quantitative explanation of the 
properties of hadrons through a solution of QCD. 

Quantum chromodynamics is the strong-interaction part of the Standard Model of 
Particle Physics and solving QCD presents a fundamental problem that is unique in the 
history of science. Never before have we been confronted by a theory whose elementary 
excitations are not those degrees-of-freedom readily accessible via experiment; i.e., 
whose elementary excitations are confined. Moreover, we have numerous reasons to 
believe that QCD generates forces which are so strong that less-than 2% of a nucleon's 



mass can be attributed to the so-called current-quark masses that appear in the QCD 
Lagrangian; viz., forces that generate mass from almost nothing, a phenomenon known 
as dynamical chiral symmetry breaking (DCSB). 

Neither confinement nor DCSB is apparent in QCD's Lagrangian and yet they play 
the dominant role in determining the observable characteristics of real-world QCD. The 
physics of hadrons is ruled by emergent phenomena, such as these, which can only be 
elucidated through the employment of nonperturbative methods in quantum field theory. 
This is both the greatest novelty and the greatest challenge within the Standard Model. 
We must find essentially new ways and means to explain precisely via mathematics the 
observable content of QCD. 

This contribution to these Proceedings provides an introduction to hadron physics and 
a review of selected recent progress in this field made using QCD's Dyson-Schwinger 
equations (DSEs). The complex of DSEs is a powerful tool, which has been employed 
with marked success to study confinement and DCSB, and their impact on hadron 
observables. This is apparent from the detailed background material that is available 
in Refs. [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1 1] and will be exemplified herein. 

2. HADRON PHYSICS - SOME KEY POINTS 

The basic problem of hadron physics is to solve QCD. This inspiring goal will only be 
achieved through a joint effort from experiment and theory. The hadron physics com- 
munity now has a range of major facilities that are accumulating data, of unprecedented 
accuracy and precision, which pose important challenges for theory. Notable amongst 
these facilities is the Thomas Jefferson National Accelerator Facility (JLab) in Newport 
News, Virginia. The interpretation and prediction of phenomena at JLab provide the 
themes for much of the material presented herein. In hadron physics, it is the feedback 
between experiment and theory that leads most rapidly to progress in understanding. 
The opportunities for researchers in hadron physics promise to grow because upgraded 
and new facilities will appear on a five-to-ten-year time-scale [12, 13, 14, 15]. 

Asymptotic coloured states have not been observed, but is it a cardinal fact that 
they cannot? No solution to QCD will be complete if it does not explain confinement. 
This means confinement in the real world, which contains quarks with light current- 
quark masses. This is distinct from the artificial universe of pure-gauge QCD without 
dynamical quarks, studies of which tend merely to focus on achieving an area law for a 
Wilson loop and hence are irrelevant to the question of light-quark confinement. 

In stepping toward an answer to the question of confinement, it will likely be nec- 
essary to map out the long-range behaviour of the interaction between light-quarks; 
namely, QCD's /3-function at infrared momenta. In this connection we note that the 
spectrum of meson and baryon excited states, and hadron elastic and transition form fac- 
tors provide unique information about the long-range interaction between light-quarks 
and, in addition, the distribution of a hadron's characterising properties - such as mass 
and momentum, linear and angular - amongst its QCD constituents. The upgraded and 
promised future facilities will provide data that should guide the charting process. How- 
ever, to make full use of that data, it will be necessary to have Poincare covariant theoret- 
ical tools that enable the reliable study of hadrons in the mass range l-2GeV. Crucially, 
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FIGURE 1. Dressed-quark mass function, M{p) in Eq. (5): solid curves - DSE results, obtained as 
explained in [17, 18], "data" - numerical simulations of lattice-QCD [19]. (NB. m = 70MeV is the 
uppermost curve and current-quark mass decreases from top to bottom.) One observes the current-quark of 
perturbative QCD evolving into a constituent-quark as its momentum becomes smaller. The constituent- 
quark mass arises from a cloud of low-momentum gluons attaching themselves to the current-quark. This 
is dynamical chiral symmetry breaking (DCSB): an essentially nonperturbative effect that generates a 
quark mass from nothing; namely, it occurs even in the chiral limit. (Figure adapted from Ref. [16].) 



on this domain both confinement and DCSB are germane. 

It is known that DCSB; namely, the generation of mass from nothing, does take place 
in QCD. It arises primarily because a dense cloud of gluons comes to clothe a low- 
momentum quark [16]. This is best seen by solving the DSE for the dressed-quark 
propagator; i.e., the gap equation, which yields the result illustrated in Fig. 1. However, 
the origin of the interaction strength at infrared momenta, which guarantees DCSB 
through the gap equation, is currently unknown. This relationship ties confinement to 
DCSB. The reality of DCSB means that the Higgs mechanism is largely irrelevant to the 
bulk of normal matter in the universe. Instead the single most important mass generating 
mechanism for light-quark hadrons is the strong interaction effect of DCSB; e.g., one 
can identify it as being responsible for 98% of a proton's mass. 

It is natural to ask whether the connection between confinement and DCSB is acci- 
dental or causal. There are models with DCSB but not confinement, however, a model 
with confinement but lacking DCSB has not yet been identified (see, e.g.. Sees. 2.1 and 
2.2 of Ref. [8]). This leads to a conjecture that DCSB is a necessary consequence of 
confinement. It is interesting that there are numerous models and theories which ex- 
hibit both confinement and DCSB, and possess an external control parameter such that 
deconfinement and chiral symmetry restoration occur simultaneously at some critical 
value of this parameter; e.g., quantum electrodynamics in three dimensions with Nf 
electrons [20, 21], and models of QCD at nonzero temperature and chemical potential 
[22, 23, 24, 25, 26]. Whether this simultaneity is a property possessed by QCD, and/or 
some broader class of theories, in response to changes in: the number of light-quark 
flavours; temperature; or chemical potential, is a longstanding question. 




The momentum-dependence of the quark mass, illustrated in Fig. 1, is an essentially 
quantum field theoretic effect, unrealisable in quantum mechanics, and a fundamental 
feature of QCD. This single curve connects the infrared and ultraviolet regimes of the 
theory, and establishes that the constituent-quark and current-quark masses are simply 
two connected points separated by a large momentum interval. The curve shows that 
QCD's dressed-quark behaves as a constituent-quark, a current-quark, or something in 
between, depending on the momentum of the probe which explores the bound-state con- 
taining the dressed-quark. It follows that calculations addressing momentum transfers 
g2 y j^2^ where M is the mass of the hadron involved, require a Poincare-covariant ap- 
proach that can veraciously realise quantum field theoretical effects [27]. Owing to the 
vector-exchange character of QCD, covariance also guarantees the existence of nonzero 
quark orbital angular momentum in a hadron's rest-frame [9, 28, 29]. 

The dressed-quark mass function has a remarkable capacity to correlate and be party 
to explaining a wide range of diverse phenomena. This brings urgency to the need to 
understand the relationship between parton properties in the light-front frame, whose 
peculiar properties simplify some theoretical analyses, and the structure of hadrons as 
measured in the rest frame or other smoothly related frames. This is a problem because, 
e.g., DCSB, an established keystone of low-energy QCD, has not been realised in the 
light- front formulation. The obstacle is the constraint := + > for massive 
quanta on the light front [30]. It is therefore impossible to make zero momentum Fock 
states that contain particles and hence the vacuum is trivial. Only the zero modes of 
light-front quantisation can dress the ground state but little progress has been made with 
understanding just how that might occur. On the other hand, it is conceivable that this 
dressing is inextricably tied with the formation and structure of Goldstone modes and not 
otherwise a measurable property of the vacuum. This conjecture is being explored [31, 
32]. In addition, parton distribution functions, which have a probability interpretation in 
the infinite momentum frame, must be calculated in order to comprehend their content: 
parametrisation is insufficient. It would be very interesting to know, e.g., how, if at all, 
the distribution functions of a Goldstone mode differ from those of other hadrons [10]. 

3. CONFINEMENT 

It is worth stating plainly that the potential between infinitely-heavy quarks measured 
in numerical simulations of quenched lattice-regularised QCD - the so-called static 
potential - is simply irrelevant to the question of confinement in the real world, in which 
light quarks are ubiquitous. In fact, it is a basic feature of QCD that light-particle creation 
and annihilation effects are essentially nonperturbative and therefore it is impossible in 
principle to compute a potential between two light quarks [33, 34]. 

A perspective on confinement was laid out in Ref. [35]. Confinement can be related to 
the analytic properties of QCD's Schwinger functions, which are often called Euclidean- 
space Green functions. For example, it can be read from the reconstruction theorem 
[36, 37] that the only Schwinger functions which can be associated with expectation 
values in the Hilbert space of observables; namely, the set of measurable expectation 
values, are those that satisfy the axiom of reflection positivity. This is an extremely 
tight constraint. It can be shown to require as a necessary condition that the Fourier 



transform of the momentum-space Schwinger function is a positive-definite function of 
its arguments, and is discussed and illustrated in Sec. 2 of Ref. [9]. 

From this standpoint the question of light-quark confinement can be translated into the 
challenge of charting the infrared behavior of QCD's universal j8 -function. (Although 
this function may depend on the scheme chosen to renormalise the theory, it is unique 
within a given scheme [38]. Of course, the behaviour of the /3 -function on the perturba- 
tive domain is well known.) This is a well-posed problem whose solution is an elemental 
goal of modem hadron physics and which can be addressed in any framework enabling 
the nonperturbative evaluation of renormalisation constants. 



4. NONPERTURBATIVE DSE TRUNCATIONS 

4.1. Existence and exact results 

Through the gap and Bethe-Salpeter equations (BSEs) the pointwise behaviour of the 
j8 -function determines the pattern of chiral symmetry breaking; e.g., the behaviour in 
Fig. 1. Moreover, the fact that these and other DSEs connect the /3-function to exper- 
imental observables entails, e.g., that comparison between computations and observa- 
tions of the hadron mass spectrum, and hadron elastic and transition form factors, can 
be used to chart the /3-function's long-range behaviour. 

In order to realise this goal, a nonperturbative symmetry-preserving DSE truncation is 
necessary. Steady quantitative progress can be made with a scheme that is systematically 
improvable [39, 40]. In fact, the mere existence of such a scheme has enabled the proof 
of exact nonperturbative results in QCD. For example, there are veracious statements 
about radially-excited and hybrid pseudoscalar mesons [41, 42]; and regarding the rj-r]' 
complex and n'^-rj-r}' mixing, with predictions of 0^1^^/ = —15° and O^^o^ = 1° [43]. In 
connection with these systems, only theoretical studies that are demonstrably consistent 
with the results proved in Refs. [41, 42, 43] can be considered seriously. 

In addition, a novel result for the pion susceptibility was recently obtained via the 
isovector-pseudoscalar vacuum polarisation [44]; viz., in the neighbourhood of the 
chiral-limit the pion susceptibility can be expressed as a sum of two independent terms: 

^5{Q = — 7^ + ^{Q+0{m) . (1) 

The first expresses the pion-pole contribution and involves the so-called vacuum-quark 
condensate. The second is identical to the vacuum chiral susceptibility, which describes 
the response of QCD's ground-state to a fluctuation in the current-quark mass. 

Equation (1) is a remarkable result, which is nonetheless readily understood. Re- 
call that in the absence of a current-quark mass, the two-flavour QCD action has a 
SUl(2)® SUi;(2) symmetry; and, moreover, that ascribing scalar-isoscalar quantum 
numbers to the QCD vacuum is a convention, contingent upon the form of the current- 
quark mass term. It follows that the massless action cannot distinguish between the 



continuum of sources specified by 



constant x J d^x q{x) e'^'^-^q{x) , 1 1 G [0, 2;r) . (2) 

Hence, the regular part of the vacuum susceptibility must be identical when measured 
as the response to any one of these sources, so that J?^eg. = ^ for all choices of d. 
This is the content of the so-called "Mexican hat" potential, which is used in building 
models for QCD. The magnitude of ^ depends on whether the chiral symmetry is 
dynamically broken, or not; and the strength of the interaction as measured with respect 
to the critical value required for DCSB [45]. When the symmetry is dynamically broken, 
then the Goldstone modes appear, by convention, in the pseudoscalar-isovector channel, 
and thus the pole contributions appear in but not in the chiral susceptibility. It is 
valid to draw an analogy with the Weinberg sum rule [46, 47]. 



4.2. Expressing DCSB in the Bethe-Salpeter kernel 

Despite the successes achieved with the systematic scheme, one anticipates that sig- 
nificant qualitative advances in understanding the essence of QCD could be made with 
symmetry-preserving kernel Ansdtze that expresses important additional nonperturba- 
tive effects, which are impossible to capture in any finite sum of contributions. Such an 
approach has recently become available [48] and is worth summarising herein. 

Consider, e.g., pseudoscalar and axial-vector mesons, which appear as poles in the 
inhomogeneous BSE for the axial- vector vertex, F^^, where f,g are flavour labels. An 
exact form of that equation is {k, q are relative momenta, P is the total momentum 
flowing into the vertex, and q± = q±P/2, etc.) 

+ J^S^^aP (k-q)^ raSf{q+) y A^^ {k, q;P), (3) 

where is a 4-point Schwinger function that is completely defined via the inverse of 
the dressed-quark propagator; namely, the kernel of the gap equation:^ 

S{p)-'=Z2{ir-p + m'''^)+Zij\^D^y{p-q)^r^S{q)^{q,p), (4) 

where indicates a Poincare invariant regularisation of the integral, with A the regular- 
isation mass-scale, D^y is the renormalised dressed-gluon propagator, is the renor- 



1 In our Euclidean metric: {7^,7v} = 25^v; vji = 7ti; 75 = 747i7273; (^fiv = Oy2)[7^,7v]; a = iti "A'; 
and Pfj timelike => < 0. More information is available, e.g., in App. A of Ref. [6]. 



malised dressed-quark-gluon vertex, and m^^ is the quark's A-dependent bare current- 
mass. The vertex and quark wave-function renormalisation constants, Zi 2(^^, A^), de- 
pend on the gauge parameter. The solution to Eq. (4) has the form 

Sip) = 1 ^(^^^^) (5) 

where the mass function, M{p^) = B{p^, C,^) /A{p^, C,^) is independent of the renormal- 
isation point, C,. The pseudoscalar vertex satisfies an analogue of Eq. (3) and has the 
general form 

iT{^{k-p) = 75 [/£5(^;^) + r-^^5(^;^) + r-^G5(fc;P) -fc7^v^^Pv//5(^;^)] ■ (6) 

In any dependable study of light-quark hadrons the solution of Eq. (3) must satisfy the 
axial-vector Ward-Takahashi; viz., 

P^T{l{k-P) + / [mf{Q +m,{Q]T{'{k;P) = S-f\k+)iys + iysS-\k^) , (7) 

which expresses chiral symmetry and its breaking pattern. The condition 

where A^^ is the analogue of in the pseudoscalar equation, is necessary and 

sufficient to ensure the Ward-Takahashi identity is satisfied [48]. 

Consider Eq. (8). Rainbow-ladder is the leading-order term in a systematic DSE 
truncation scheme [39, 40]. It corresponds to r{ = jy, in which case Eq. (8) is solved by 
Ag^^ = = Aj|. This is the solution that indeed provides the rainbow-ladder forms 
of Eq. (3). Such consistency will be apparent in any valid systematic term-by-term 
improvement of the rainbow-ladder truncation. 

However, Eq. (8) is far more than just a device for checking a truncation's consistency. 
For, just as the vector Ward-Takahashi identity has long been used to build Ansdtze for 
the dressed-quark-photon vertex [1, 49, 50], Eq. (8) provides a tool for constructing a 
symmetry-preserving kernel of the BSE that is matched to any reasonable Ansatz for 
the dressed-quark-gluon vertex which appears in the gap equation. With this powerful 
capacity, Eq. (8) achieves a goal that has been sought ever since the Bethe-Salpeter 
equation was introduced [51]. As we shall indicate, the symmetry -preserving kernel 
it can provide promises to enable the first reliable Poincare invariant calculation of the 
spectrum of mesons with masses larger than 1 GeV. 

One can illustrate the utility of Eq. (8) through an application to ground state pseu- 
doscalar and scalar mesons composed of equal-mass u- and J-quarks. To this end, sup- 
pose that in Eq. (4) one employs an Ansatz for the quark-gluon vertex which satisfies 



P^iTf^{k+,k^)=^{P^) S-Mk+)-S-r\k^ 



(9) 



with flavour-independent. (NB. While the true quark-gluon vertex does not satisfy 
this identity, owing to the form of the Slavnov-Taylor identity which it does satisfy, it is 



plausible that a solution of Eq. (9) can provide a reasonable pointwise approximation to 
the true vertex.) Given Eq. (9), then Eq. (8) entails {I = q — k) 



ilpK{l{k,q-P) = ^[If \r{\q-P)-T{\k-P) 



(10) 



with an analogous equation for P^l piA^^^p{k, q; P) . This identity can be solved to obtain 



A-^^l{k,q;P):=^{{k-q)^)r5A^l^'{k,q;P), (11) 



v^'ith, using Eq. (6), 



A-jf {k,q;P) = 2£p [iAe, {q,k-P) + y- PAp, iq,k;P)] + Yp k; P) 

+Upy-lAGM^k-P) + [yp,y-P]I.HM^k-P)+2lp[y-l,y-P]AHM^k-P)An) 

where £ - (q + k)/2, Zct>{q,k;P) = [^{q;P) +^{k;P)]/2 and A<^{q,k\P) = [^{q;P) - 
c^{k;P)]/[q'^-k^]. 

Now, given any Ansatz for the quark-gluon vertex that satisfies Eq. (9), then the pseu- 
doscalar analogue of Eq. (3) and Eqs.(4), (11), (12) provide a symmetry-preserving 
closed system whose solution predicts the properties of pseudoscalar mesons. The rele- 
vant scalar meson equations are readily derived. (NB. We are aware of the role played 
by resonant contributions to the kernel in the scalar channel [52] but they are not perti- 
nent to this discussion.) With these systems one can anticipate, elucidate and understand 
the influence on hadron properties of the rich nonperturbative structure expected of the 
fully-dressed quark-gluon vertex in QCD: in particular, that of the dynamically gener- 
ated dressed-quark mass function, whose impact is quashed at any finite order in the 
truncation scheme of Ref . [40] . 

To proceed one need only specify the gap equation's kernel because, as noted above, 
the BSEs are completely defined therefrom. To complete the illustration [48] a simplified 
form of the effective interaction in Ref. [53] was employed and two vertex An^aYze were 
compared; viz., the bare vertex = 7^, which defines the rainbow-ladder truncation 
of the DSEs and omits vertex dressing; and the Ball-Chiu (BC) vertex [49] which 
nonperturbatively incorporates vertex dressing associated with DCSB: 

/r^(^,fc) = £a, (^2^^2)7^ +2£^ [iy-lAAs{q^,k^)+ABs{q^.k^)\. (13) 

A particular novelty of the study is that one can calculate the current-quark-mass- 
dependence of meson masses using a symmetry-preserving DSE truncation whose dia- 
grammatic content is unknown. That dependence is depicted in Fig. 2 and compared with 
the rainbow-ladder result. The m-dependence of the pseudoscalar meson's mass provides 
numerical confirmation of the algebraic fact that the axial- vector Ward-Takahashi iden- 
tity is preserved by both the rainbow-ladder truncation and the BC-consistent Ansatz for 
the Bethe-Salpeter kernel. The figure also shows that the axial-vector Ward-Takahashi 
identity and DCSB conspire to shield the pion's mass from material variation in response 
to dressing the quark-gluon vertex [8, 54]. 
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FIGURE 2. Dependence of pseudoscalar (left panel) and scalar (right) meson masses on the current- 
quark mass, m. The Ball-Chiu vertex (BC) result is compared with the rainbow-ladder (RL) result. (Figure 
adapted from Ref. [48].) 



As noted in Ref. [48], a rainbow-ladder kernel with realistic interaction strength yields 



2M(0) 



(0.3±0.i; 



(14) 



RL 



which can be contrasted with the value obtained using the BC-consistent Bethe-Salpeter 
kernel; viz., 

(15) 



eBC<o.l. 



Plainly, significant additional repulsion is present in the BC-consistent truncation of the 
scalar BSE. 

Scalar mesons are commonly identified as ^Pq states. This assignment expresses 
a constituent-quark-model perspective, from which a 7^^ = 0++ fermion-antifermion 
bound-state must have the constituents' spins aligned and one unit of constituent orbital 
angular momentum; and hence a scalar is a spin and orbital excitation of a pseudoscalar 
meson. Of course, no constituent-quark-model can be connected systematically with 
QCD. Nevertheless, as we observed above, the presence of orbital angular momentum 
in a hadron's rest frame is a necessary consequence of Poincare covariance and the 
vector-boson-exchange character of QCD [9, 28, 29], so there is a realisation in QCD of 
the quark-model anticipation. 

Extant studies of realistic corrections to the rainbow-ladder truncation show that 
they reduce hyperfine splitting [54]. Hence, with the comparison between Eqs.(14) 
and (15) one has a clear indication that in a Poincare covariant treatment the BC- 
consistent truncation magnifies spin-orbit splitting, an effect which can be attributed to 
the influence of the quark's dynamically-enhanced scalar self-energy [9] in the Bethe- 
Salpeter kernel. 
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FIGURE 3. Left panel - solution of vector-vertex Bethe-Salpeter equation; and right panel - analogue 
for the axial-vector vertex. In both panels, short-dashed curve - rainbow-ladder truncation; long-dashed 
curve - Ball-Chiu-consistent truncation; and solid curve - improved BC-vertex, with quark anomalous 
chromomagnetic moment, Eqs. (16), (17). From the zero in the curves, one reads the mass-squared of the 
lightest meson supported by the assumed kernel in each channel. 

5. QUARK ANOMALOUS CHROMOMAGNETIC MOMENT 

It was conjectured in Ref. [48] that the full realisation of DCSB in the Bethe-Salpeter 
kernel will have a material impact on mesons with mass greater than 1 GeV. Moreover, 
that it can overcome a longstanding failure of theoretical hadron physics. Namely, 
no extant hadron spectrum calculation is believable because all symmetry preserving 
studies produce a splitting between vector and axial-vector mesons that is far too small: 
just one-quarter of the experimental value (see, e.g., Refs. [55, 56, 57]). 

In this connection, preliminary results^ are now available [58] and they are sum- 
marised in Fig. 3. The curves are obtained by solving the inhomogeneous vector- and 
axial-vector-vertex BSEs at spacelike momenta, and fitting a Pade approximant to ex- 
trapolate the result to timelike momenta. Each vertex has a dominant Dirac covariant: 
YfiFy{k;P), y^Fl{k;P), respectively; and the curves in Fig. 3 depict the -dependence 
of \/Fyj^{k = 0;P^) =: Gv,a{P^)- The location of the first zero in such a curve yields 
the mass-squared of the lightest bound-state produced by the assumed interaction in the 
associated channel [59]. The masses thus determined are listed in Table 1. 

To illustrate the longstanding difficulty, using the interaction described in Ref. [48], 
the masses were first calculated in rainbow-ladder truncation; viz., leading-order in the 
systematic and symmetry-preserving truncation scheme of Ref. [40] . As anticipated, 
while the p-meson mass is acceptable, the ai-mass is far too small. 

As we saw above, the procedure introduced in Ref. [48] enables meson masses to be 
calculated using a symmetry-preserving DSE truncation whose diagrammatic content 
is unknown. One can therefore elucidate the effect of an essentially nonperturbative 
DCSB component in dressed-quark gluon vertex on the p-ai complex; in this case, the 
Ag term in Eq. (13) - recall the enormous impact of this term in the scalar channel. 



Preliminary in the sense that not all Dirac covariants are employed in solving for the vertices and internal 
consistency checks for the pseudoscalar and scalar channels are still underway. 



TABLE 1. Axial-vector and vector meson masses calculated in three truncations of the coupled 
gap and Bethe-Salpeter equations. The last column was obtained using the standard Ball-Chiu An^afz 
augmented by the quark anomalous chromomagnetic moment in Eqs. (16), (17). 
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Fig. 2. The result is shown in Fig. 3: the DCSB Ag-term boosts the ai mass, which is a 
positive outcome, but it simultaneously boosts the p mass, such that the mass-splitting 
is practically unchanged from the rainbow-ladder result - see Table 1 . Was Ref. [48] 
too optimistic in expressing a hope that the scheme introduced could provide the first 
realistic meson spectrum that also encompassed states with mass greater-than 1 GeV? 

Before answering, let us return to a consideration of chirally symmetric QCD. That 
theory exhibits helicity conservation so that, perturbatively, the quark-gluon vertex can- 
not have a term with the helicity-flipping characteristics of A5. There is another feature 
of massless fermions in gauge field theories; namely, they cannot posses an anoma- 
lous chromo/electro-magnetic moment because the term that describes it couples left- 
and right-handed fermions. However, if the theory's chiral symmetry is strongly bro- 
ken dynamically, why shouldn't the fermions posses a large anomalous chromo/electro- 
magnetic moment? Such an effect is expressed in the quark-gluon-vertex via a term 

r'^'^{q,k)^G^y{q-k)yT5{q.k) (16) 

where, owing to DCSB, the natural strength is represented by theAnsatz 

T5{q,k)=AB{q\k^). (17) 

NB. Based on the models in Refs. [27, 60], 2Mi7Z(M|)As(M|,M|) ~ where Me is 
the Euclidean constituent-quark mass, defined in Ref. [53]. 

Using the procedure introduced in Ref. [48], the vector and axial- vector vertex equa- 
tions can be solved using the dressed-quark-gluon vertex obtained as the sum of 
Eqs. (13) and (16). The effect is remarkable and plain in Fig. 3: the anomalous chro- 
momagnetic moment leads to additional repulsion in the ai channel but significant at- 
traction in the p channel such that, for the first time, a realistic result is simultaneously 
obtained for the masses in both these channels, and hence the ai-p mass- splitting - see 
Table 1 . Furthermore, the origin of the splitting, in an interference between Ag in Eq. (13) 
and T5 in Eqs. (16), (17) is intuitively appealing. In the chiral limit the mass-splitting be- 
tween parity partners should owe solely to DCSB and here that is seen explicitly: in the 
absence of DCSB, Ag = = T5. The rainbow-ladder result is also understood: this trun- 
cation fails to adequately express DCSB in the Bethe-Salpeter kernel and hence cannot 
realistically split parity partners. 

Notably, Table 1 is a "first-guess" result; i.e., there was no tuning of the strength in 
Eq. (17), so how reliable can it be? This question amounts to deciding whether a realistic 



size is assumed for a light-quark's anomalous chromomagnetic moment. Fortunately, 
an analysis is available of results for the dressed-quark-gluon vertex obtained through 
numerical simulations of quenched-QCD [61]. This study shows that T5 is dynamically 
two orders -of -magnitude larger than the one-loop perturbative result and, indeed, is of 
the same magnitude and possesses the same domain of significant support as lS.B{q^, k^), 
precisely in accordance with the assumption we have made. 

At this point it is natural to consider whether DCSB in QCD can also generate a 
large quark anomalous electromagnetic moment term in the quark-photon vertex. In 
perturbation theory, of course, since it doesn't express DCSB, the quark's anomalous 
electromagnetic moment is small [62]. One obtains the same answer in the rainbow- 
ladder truncation; e.g., the Fg and Fg terms in Ref. [63], which combine to form T5 in our 
notation, contribute less-than 1% to the pion's electromagnetic form factor. However, 
as we've already seen, this truncation doesn't adequately incorporate DCSB into the 
Bethe-Salpeter kernel. At present there is no concrete information available that can be 
used to answer the question but that is about to change. Upon completion of the study of 
the ai-p complex, we will have determined the magnitude of the analogue of T5 in the 
colour-singlet vector vertex. It could plausibly be large. 

We reiterate that the large quark anomalous chromomagnetic moment owes to DCSB, 
which itself is probably a consequence of confinement; and real-world electrons and 
muons do not couple directly to interactions associated with, or degrees-of-freedom 
affected by, DCSB. There is thus little reason to anticipate measurable contributions 
to the electron and muon anomalous electromagnetic moments. Nevertheless, given the 
magnitude of the muon "^^ — 2 anomaly" and its assumed importance as an harbinger 
of physics beyond the Standard Model [64], it might be worthwhile to compute a 
quantitative estimate of the contribution to = (g^ — 2)/2 from the quark's DCSB- 
induced anomalous chromomagnetic, and possibly electromagnetic, moments. Such 
contributions will appear in the hadronic component of the photon polarisation tensor. 

6. PION ELECTROMAGNETIC FORM FACTOR 

In charting the long-range interaction between light-quarks via the feedback between 
experiment and theory, hadron elastic and transition form factors can provide unique in- 
formation, beyond that obtained through studies of the hadron spectrum. This is demon- 
strated very clearly by an analysis of the electromagnetic pion form factor, F^^{Q^), 
because the pion has a unique place in the Standard Model. It is a bound-state of a 
dressed-quark and -antiquark, and also that almost-massless collective excitation which 
is the Goldstone mode arising from the dynamical breaking of chiral symmetry. This 
dichotomy can only be understood by merging the study of many-body aspects of the 
QCD vacuum with the symmetry-preserving analysis of two-body bound-states [65]. 
Furthermore, the possibility that this dichotomous nature could have wide-ranging ef- 
fects on pion properties has made the empirical investigation of these properties highly 
desirable, despite the difficulty in preparing a system that can act as a pion target and the 
concomitant complexities in the interpretation of the experiments; e.g., [66, 67, 68, 69]. 

The merit of using F^^{Q^) to elucidate the potential of an interplay between exper- 
iment and nonperturbative theory, as a means of constraining the long-range behaviour 



of QCD's /3-function, is amplified by the existence of a prediction [70, 71, 72] that 
Q^Pit{Q^) ~ constant for ^ m| in a theory whose interaction is mediated by mass- 
less vector-bosons. The verification of this prediction is a strong motivation for modern 
experiment [66, 67, 68], which can also be viewed as an attempt to constrain and map 
experimentally the pointwise behaviour of the exchange interaction that binds the pion. 

Poincare covariance entails that the Bethe-Salpeter amplitude for an isovector pseu- 
doscalar bound-state of a dressed-quark and -antiquark takes the form 

(18) 

where {t', j = 1,2,3} are the Pauli matrices. This amplitude is determined from a 
homogeneous BSE; exemplified, e.g., by Eq. (3) with the driving term, 22757^, omitted. 

We can now explicate one of the most important consequences of DCSB. Namely, a 
set of four Goldberger-Treiman relations, exact in chiral-limit QCD [65]: 

jiE^{k-0)=Bo{k^), FR{k;0)+2f^F^{k;0) =Ao{k^), 

GR{k;0)+2jiGn{k;0)=A'o{k^), HR{k;0) +2f^H^{k;0) =0. ^ 

Here AQ{k^), BQ{k^) describe the solution of the chiral limit gap equation; is the 
pion's chiral-limit leptonic decay constant; and one has used the following features of 
the chiral-limit axial- vector vertex: 

+y75 [rnFn{k-P) + Y^Y-kGn{k;P)-o^ykyHn(k;P)] +rl'^{k;P), (20) 

where Fr,Gr, Hr, T'^^ {k; P) are regular for P^c^O and P^T"^^ {k; P) ~ 0(/'2) for P^ ~ 0. 

Equations (19) are fascinating and have far-reaching consequences. The first identity 
states that, owing to DCSB, the dominant piece of the pseudoscalar meson Bethe- 
Salpeter amplitude is completely determined by the scalar piece of the dressed-quark 
self-energy; i.e., by solving the quark one-body problem, one obtains simultaneously 
the solution for an important part of the pseudoscalar meson two-body problem. The 
next two demonstrate that a pseudoscalar meson necessarily contains components of 
pseudovector origin. These terms alter the asymptotic form of F^"^{Q^) by a factor of 
cf. the result obtained in their absence [73]. One may consider Eqs. (19) as an expression 
of pointwise consequences of Goldstone's theorem, which is therefore seen to have an 
impact far broader than that described in textbooks. 

QCD-based DSE calculations of F^^{Q^) exist [73, 74], the most systematic of which 
[74] predicted the measured form factor [66]. Germane to our discourse, however, is an 
elucidation of the sensitivity of F^^{Q^) to the pointwise behaviour of the interaction 
between quarks. We therefore recapitulate on Ref. [75], which explored how predictions 
for pion properties change if quarks interact not via massless vector-boson exchange but 
instead through a contact interaction; viz.. 



(21) 



TABLE 2. Results obtained with (in GeV) m = 0, mc = 0.1 1 , Air = 1 h\x = 0.24, 
Auv = l/fuv = 0.823. They are commensurate with those from QCD-based DSE 
studies [53]. Dimensioned quantities are Usted in GeV or fm, as appropriate, and 

Fji Fr M K fn\f^^Q ''ii ''iilp^^Q 

0.23 4.28 0.69 0.68 0.40 0.22 0.094 0.11 0.29 0.41 



where mq is a gluon mass-scale (such a scale is generated dynamically in QCD, with a 
value ~ 0.5 GeV [76]), and proceeded by embedding this interaction in a rainbow-ladder 
truncation of the DSEs. 

In this case, using a confining regularisation scheme [77], the gap equation, which de- 
termines this interaction's momentum-independent dressed-quark mass, can be written 



M = m + 



M 



(22) 



where m is the current-quark mass and jM = r( — 1,M T„y) — r( — 1,M Tj^), with 
r{a,y) being the incomplete gamma-function. Results are presented in Table 2. 

With a symmetry -preserving regularisation of the interaction in Eq. (21), the Bethe- 
Salpeter amplitude cannot depend on relative momentum. Hence Eq. (18) becomes 



1 



r^iP) = 75 

and the explicit form of the model's ladder BSE is 



(23) 
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where, with m = = P^, anticipating the Goldstone character of the pion, 

2J(rEE = 1^i{M\xlxl 



J^pp 



EE ■ 



(25) 



and ^i(z) = —z'io"{z), where we have suppressed the dependence on Tir uv The solution 
of Eq. (24) gives the pion's chiral-limit Bethe-Salpeter amplitude: 



eI = QM1,fI =0.160, 



(26) 



written with unit normalisation. 

However, this is not the physical convention. The canonical normalisation procedure 
ensures unit residue for the pion bound-state contribution to the quark-antiquark scatter- 
ing matrix, a property of T^^{P) = -y r]j.(P), where 



"^'^ ri{-P)4rSiq+P)^iiP)s{q). 



{2n 



dPu 



(27) 



In the chiral limit, 



^0 = ^^^x{m\xItI)e],[eI-2F},]. (28) 

The pion's leptonic decay constant is obtained from the canonically normalised ampli- 
tude and in the chiral limit 

If one has preserved Eq. (7), then, for m = in the neighbourhood of = 0, the 
solution of the axial- vector BSE has the form - cf. Eq. (20): 

and the following subset of Eqs. (19) will hold: 

f^E'^=M,2^+FR = l. (31) 

That they do can be verified from Table 2, which also shows that Fti{P), necessarily 
nonzero in a vector exchange theory, irrespective of the pointwise behaviour of the 
interaction, has a measurable impact on the value of fn- 

Based upon these results, one can proceed to compute the electromagnetic pion 
form factor in the generalised impulse approximation [75]; i.e., at leading-order in a 
symmetry-preserving DSE truncation scheme [73, 74, 78]. Namely, for an incoming 
pion with momentum pi =K — Q/2, which absorbs a photon with space-like momentum 
Q, so that the outgoing pion has momentum p2 = K-\-Q/2, 



2K,Fr{Q')=2N,.J^^tr^ 



nii-P2)sit+p2)irfiSit+pi) s{t) 



(32) 



i2ny 

The form factor is expressible as a sum; viz., 

^r(e') = F^^EE{Q')+F^':'EFiQ')+FrFFiQ'), (33) 
= E^,^TMQ')+ElF^T-,{Q^)+F^%^,iQ^), (34) 

wherein each function has a simple algebraic form in this model. 

In the left panel of Fig. 4 we present F^^{Q^) and the three separate contributions 
defined in Eq. (34). We highlight two features. First, F^^p{Q^ = 0) contributes roughly 
one-third of the pion's unit charge. This could have been anticipated from Eq. (28). 
Second, and perhaps more dramatic: the interaction in Eq. (21) generates 

F^'^iQ- ^ oo) = constant. (35) 

Both results originate in the nonzero value of F^{P), which is a straightforward conse- 
quence of the symmetry -preserving treatment of a vector exchange theory [65]. Equa- 
tion (35) should not come as a surprise: with a symmetry-preserving regularisation of the 
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FIGURE 4. Left panel . F™{Q^) and the separate contributions introduced in Eq. (34). F™{Q^ = 0) = 
1, without fine-tuning, because a symmetry-preserving regularisation of the interaction in Eq. (21) was 
implemented. Right panel. Solid curve: Q^Fji{Q^) obtained with Eq. (21). Dashed curve: DSE prediction 
[74], which employed a momentum-dependent renormalisation-group-improved gluon exchange inter- 
action. For > O.lTGeV^ w M^, marked by the vertical dotted line, the contact interaction result for 
F^''\Q^) differs from that in Ref. [74] by more than 20%. The data are from Refs. [66, 67, 68]. (Figure 
adapted from Ref. [75].) 



interaction in Eq. (21), the pion's Bethe-Salpeter amplitude cannot depend on the con- 
stituent's relative momentum. This is characteristic of a pointlike particle, which must 
have a hard form factor. 

A striking outcome, evident in Fig. 4: F^^^ is only a good approximation to the net 

pion form factor for < M^. F^^^ and F^^p evolve with equal rapidity - there is no 
reason for this to be otherwise, as they are determined by the same mass-scales - but a 
nonzero constant comes quickly to dominate over a form factor that falls swiftly to zero. 

In the right panel of Fig. 4 we compare the form factor computed from Eq. (21) with 
contemporary experimental data [66, 67, 68] and a QCD-based DSE prediction [74]. 
Both the QCD-based result and that obtained from the momentum-independent inter- 
action yield the same values for the pion's static properties. However, for > the 
form factor computed using ~ 1/^^ vector boson exchange is immediately distinguish- 
able empirically from that produced by a momentum-independent interaction. Indeed, 
the figure shows that for F^^, existing experiments can already distinguish between dif- 
ferent possibilities for the quark-quark interaction. 

This summary illustrates that when a momentum-independent vector-exchange inter- 
action is regularised in a symmetry-preserving manner, the results are directly compara- 
ble with experiment, computations based on well-defined and systematically-improvable 
truncations of QCD's DSEs [74], and perturbative QCD. In this context it will now be 
apparent that a contact interaction, whilst capable of describing pion static properties 
well. Table 2, generates a form factor whose evolution with deviates markedly from 
experiment for > 0.17GeV^ ~ and produces asymptotic power-law behaviour, 
Eq. (35), in serious conflict with perturbative- QCD [70, 71, 72]. 

In closing this section we note that the contact interaction produces a momentum- 
independent dressed-quark mass function, in contrast to QCD-based DSE studies [9, 18] 



and lattice-QCD [19]. This is fundamentally the origin of the marked discrepancy be- 
tween the form factor it produces and extant experiment. Hence Ref. [75] highlights that 
form factor observables, measured at an upgraded Jefferson laboratory, e.g., are capa- 
ble of mapping the running of the dressed-quark mass function. Efforts are underway to 
establish the signals of the running mass in baryon elastic and transition form factors. 

7. PION AND KAON VALENCE-QUARK DISTRIBUTIONS 

The past forty years has seen a tremendous effort to deduce the parton distribution func- 
tions of the most accessible hadrons - the proton, neutron and pion. There are many rea- 
sons for this long sustained and thriving interest [10] but in large part it is motivated by 
the suspected process-independence of the usual parton distribution functions and hence 
an ability to unify many hadronic processes through their computation. In connection 
with uncovering the essence of the strong interaction, the behaviour of the valence-quark 
distribution functions at large Bjorken-x is most relevant. Importantly, in the infinite mo- 
mentum frame, Bjorken-x measures the fraction of a hadron's four-momentum carried 
by the struck parton and, e.g., the valence-quark distribution function, qy{x), measures 
the number-density of valence-quarks with momentum-fraction x. 

Owing to the dichotomous nature of Goldstone bosons, understanding the valence- 
quark distribution functions in the pion and kaon is of great importance. Moreover, 
given the large value of the ratio of s-io-u current-quark masses, a comparison between 
the pion and kaon structure functions offers the chance to chart effects of explicit chiral 
symmetry breaking on the structure of would-be Goldstone modes. There is also the 
prediction [79, 80] that a theory in which the quarks interact via l/k^ vector-boson 
exchange will produce valence-quark distribution functions for which 

^vWoc(l-x)2+r,x>0.85, (36) 

where 7 > is an anomalous dimension that grows with increasing momentum transfer. 
(See Sec.VI.B.3 of Ref. [10] for a detailed discussion.) 

Experimental knowledge of the parton structure of the pion and kaon arises primarily 
from pionic or kaonic Drell-Yan scattering from nucleons in heavy nuclei [69, 81, 82]. 
Theoretically, given that DCSB plays a crucial role in connection with pseudoscalar 
mesons, one must employ an approach that realistically expresses this phenomenon. 
The DSEs therefore provide a natural framework: a study of the pion exists [83] and one 
of the kaon is underway [84]. 

One can illustrate the results anticipated from the latter study through an internally 
consistent calculation based upon the QCD-improvement of a simple model used al- 
ready for pion and kaon distribution functions [85]. In its original form, as a version 
of the Nambu-Jona-Lasinio model with a hard-cutoff, A = 0.9GeV, and momentum- 
independent dressed-quark masses, M„ = 0.35 GeV and Ms = 0.53 GeV, the model yields 
valence-quark distribution functions that behave as ( 1 — x) at large x, in conflict with per- 
turbative QCD and nonperturbative DSE studies. However, this can be remedied through 
the inclusion of a pion Bethe-Salpeter amplitude whose dependence on relative momen- 



TABLE 3. Leading moments of our valence M-quark distributions; 
(x")g := /q dxx"u^{x; g), with Qq = 0.57 GeV and g 5 GeV. For com- 
parison, the first numerical column shows results from Ref. [83]. 
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turn is QCD-like; viz., ^nik^) °^ l/k^ when the magnitude of is large. To be explicit, 

Ik^, {k-Mu{k) 
where Kjr = -x{MI -m\{ \ - x)) /{ \ - x), with Kir - K"uv = 



<^);r ( JC) = ^^P('^/^p ^ . M„ ( k) = m„ + (M„ - m„) <I>^{k) , (38) 

which is an Ansatz for the momentum-dependent dressed-quark mass, Mi,{k^), with 
= 5.5 Me V, in a form motivated by the first of Eqs. (19); and the value of the pion- 
quark-antiquark coupling constant, g^^^^^ ensures 

1 

dxu'^{x) = \. (39) 







A value of A;;; = 0.6GeV provides a satisfactory least-squares fit to {u^{x),x > 0.4} in 
Ref. [83]. (NB. The hard cutoff in Eq. (37), mimics the effect of a second Bethe-Salpeter 
amplitude, which would otherwise appear.) 
Repeating the analysis for the kaon, we obtain 



"fW = -^k/ dK ,;,:^2^2 -^k{.^)^ (40) 



where Kir = ^^(^.? ^ '"i ( 1 ~ ^) ) / ( 1 ~ ^) ' with Kir — '^uv = ; 

= ^^P('^/^i)-\ M,(K) =m,-F(M,-m,)<^)^(K), (41) 

with rus = 0.135 GeV and = 1.1 A;i-, the latter value motivated by the BSE results in 
Ref. [53]; and the value of gKsu ensures the kaon analogue of Eq. (39). 

In Table 3 we report low moments of the computed valence-quark distribution func- 
tions. The near numerical agreement between entries in the first two columns is an il- 
lustration of the similarity between our model distribution and that by which it was con- 
strained [83]. However, it should be borne in mind that these low moments are not very 
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FIGURE 5. Left panel. Valence M-quark distribution functions, computed as described in the text: solid 
curve, Mv(x); and dashed curve, u^{x). Applying leading-order QCD evolution from Qq — 0.32GeV^ 
to = 25GeV^, explained in Sec.II.D of Ref. [10], one obtains the other two curves from these 
starting distributions: dashed curve, u^{x;25); and dotted curve, uf{x,25). Right panel. Model ratio 
Uy /uy evaluated at Qq = 0.32GeV^, Dashed curve, and = 25GeV^, solid curve, /m" = 0.13. 
Under the right conditions, u^ /u^ should equal the ratio of kaon-to-pion Drell-Yan cross-sections, and 
we reproduce that obtained from a sample of dimuon events with invariant mass 4.1 < M < 8 .5 GeV [81]. 

sensitive to the pointwise behaviour of distributions near x = I - see, e.g., Sees. VI.B.2 
and VLB. 3 of Ref. [10]. 

In Fig. 5 we depict our computed distributions themselves and relevant ratios. As- 
pects of the curves are model-independent. For example, owing to its larger mass, one 
anticipates that the 5-quark should carry more of the charged-kaon's momentum than the 
M-quark. This explains, in a manner which is transparent in this model, why the support 
of xu^ (x) is shifted to lower-x than that xu^ (x) . QCD evolution is an area-conserving op- 
eration on the distribution function, which shifts support from large-x to small-x. Thus, 
while both Uy'^{x; 2o) °^ ( 1 - for x ~ 1 , 



These observations explain the qualitative behaviour of the evolved distributions. Con- 
cerning the ratio, as a consequence of the form of the evolution equations, the value of 
the ratio at jc = 1 is invariant under evolution. At the other extreme, the value of the ratio 
at jc = approaches one under evolution owing to the increasingly large population of 
sea-quarks produced thereby. 



While a symmetry-preserving description of mesons is essential, it is only part of the 
physics that nonperturbative QCD must describe since Nature also presents baryons, 
light-quarks in three-particle composites. An explanation of the spectrum of baryons 
and the nature of interactions between them is basic to understanding the Standard 



(42) 



8. BARYON PROPERTIES 



Model. The present and planned experimental programmes at JLab, and other facilities 
worldwide, are critical elements in this effort. 

No approach to QCD is comprehensive if it cannot provide a unified explanation of 
both mesons and baryons. We have explained that DCSB is a keystone of the Standard 
Model, which is evident in the momentum-dependence of the dressed-quark mass func- 
tion - Fig. 1: it is just as important to baryons as it is to mesons. The DSEs furnish the 
only extant framework that can simultaneously connect both meson and baryon observ- 
ables with this basic feature of QCD, having provided, e.g., a direct correlation of meson 
and baryon properties via a single interaction kernel, which preserves QCD's one-loop 
renormalisation group behaviour and can systematically be improved [86, 87]. 

In quantum field theory a baryon appears as a pole in a six-point quark Green function. 
The residue is proportional to the baryon's Faddeev amplitude, which is obtained from a 
Poincare covariant Faddeev equation that sums all possible exchanges and interactions 
that can take place between three dressed-quarks. A tractable Faddeev equation for 
baryons [88] is founded on the observation that an interaction which describes colour- 
singlet mesons also generates nonpointlike quark-quark (diquark) correlations in the 
colour-3 (antitriplet) channel [89]. The lightest diquark correlations appear in the = 
0+, 1+ channels [90, 91] and hence only they need be retained in approximating the 
quark-quark scattering matrix that appears as part of the Faddeev equation [27, 87]. 

Diquarks do not appear in the strong interaction spectrum [40, 54] but the attraction 
between quarks in this channel justifies a picture of baryons in which two quarks within 
a baryon are always correlated as a colour-3 diquark pseudoparticle, and binding is ef- 
fected by the iterated exchange of roles between the bystander and diquark-participant 
quarks. Here it is important to emphasise strongly that QCD supports nonpointlike di- 
quark correlations [92] . Hence models that employ pointlike diquark degrees of freedom 
have little connection with QCD. 

Numerous properties of the nucleon have been computed using the Faddeev equation 
just described. Herein we draw only two examples from comprehensive analyses of 
nucleon electromagnetic form factors [27, 34, 93, 94]. To introduce the results, we note 
that the nucleon's electromagnetic current is 

MP',P) = ieu{P') (^rnFi{Q^) + :^o^yQyF2iQ^)^u{P), (43) 

where P (P') is the momentum of the incoming (outgoing) nucleon, Q = P' — P is the 
momentum transfer, and Fi and F2 are, respectively, the Dirac and Pauli form factors, 
from which one obtains the nucleon's electric and magnetic (Sachs) form factors 

Ge{Q^) =Fi(e2)__^F2(e2), Gm{Q^) =Fi{Q^)+F2{Q^). (44) 

In the left panel of Fig. 6 we depict a ratio of flavour- separated contributions to the 
neutron's Dirac form factor. The predicted 2^-dependence owes to the presence of axial- 
vector diquark correlations in the nucleon. It has been found [27] that the neutron's 
singly-represented w-quark is more likely to be struck in association with an axial- 
vector diquark correlation than with a scalar, and form factor contributions involving 
an axial-vector diquark are soft. On the other hand, the doubly-represented J-quark is 
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FIGURE 6. Left panel. Computed ratio of flavour-separated contributions to neutron's Dirac form 
factor (M =nucleon mass) solid curve. Data, reconstructed from selected neutron electric form factor 
data using parametrisations of Ref. [95] when necessary [96]: Ref. [97], up-triangles; Ref. [98], diamonds; 
Ref. [99], circles; and Ref. [100], down-triangles. Right panel. Neutron magnetic form factor divided by 
dipole fit: long-dashed curve, DSE prediction [27]; solid curve, [95]; opeti circles, recent CLAS results 
[101] (red) and pre-2008 data (green); and black squares, anticipated Q^-coverage and errors for an 
experiment at the upgraded JLab [102] - the magnitude of these points is naturally arbitrary at this time 
but they were placed to follow a recent fit to now-extant experimental data {short-dashed curve, [103]). 



predominantly linked with harder scalar-diquark contributions. NB. When isospin is a 
good symmetry, f"'"/f^^''^ = F^''^ /F^'" . 

In the right panel of Fig. 6 we depict computed using the Faddeev equation [27]. 
Given that it is common practice to compare nucleon form factors with an empirical 
dipole, the DSE result is presented in a similar way. Namely, the function 1/(1 + 
was fitted to the DSE result on 2 < < 14. This domain excludes 

the region whereupon pion cloud effects are significant and maximises coverage of 
the domain on which the quark-core calculation is most reliable. The fit produced 
niD = l.OSMyv cf. the empirical value, = 0.90 Mf^f, which is just 14% smaller. 

The DSE predictions in Fig. 6 depend sensitively on the momentum-dependence of 
the dressed-quark mass function, which feeds into properties such as the diquarks' mass 
and size. Since in the computation of elastic and transition form factors, the probe's 
input momentum Q is principally shared equally amongst the dressed-quarks, then each 
absorbs a momentum fraction Q/3. Thus, to scan the behaviour of the mass function on 
the domain e [0.2, 1.2] GeV^, upon which its character changes from nonperturbative 
to perturbative, one requires E [2, 15]GeV^. The figure illustrates that comparison 
and feedback between DSE results and forthcoming precision data on nucleon form 
factors can serve as a means by which to empirically chart the momentum evolution 
of the dressed-quark mass function, and therefrom the infrared behavior of QCD's 
j8 -function. In particular, it should enable the unambiguous location of the transition 
boundary between the constituent- and current-quark domains that is signalled by the 
sharp drop apparent in Fig. 1 and which can likely be related to an infrared inflexion 
point in QCD's running coupling, whose properties are determined by the /3 -function. 



9. EPILOGUE 



Dynamical chiral symmetry breaking (DCSB) is a fact in QCD. It is manifest in dressed- 
propagators and vertices, and, amongst other things, it is responsible for: the transforma- 
tion of the light current-quarks in QCD's Lagrangian into heavy constituent-like quark's, 
in terms of which order was first brought to the hadron spectrum; the unnaturally small 
values of the masses of light-quark pseudoscalar mesons; the unnaturally strong cou- 
pling of pseudoscalar mesons to light-quarks - g^qq ~ 4.3; and the unnaturally strong 
coupling of pseudoscalar mesons to the lightest baryons - ^;i.yvA' ~ 12.8 ~ ^gnqq- 

Herein we have illustrated the dramatic impact that DCSB has upon observables: the 
spectrum. Sees. 4 and 5; hadron form factors. Sees. 6 and 8; and parton distribution func- 
tions. Sec. 7. A "smoking gun" for DCSB is the behaviour of the dressed-quark mass 
function. The momentum dependence manifest in Fig. 1 is an essentially quantum field 
theoretical effect. Exposing and elucidating its effects therefore requires a nonpertur- 
bative and symmetry-preserving approach, where the latter means preserving Poincare 
covariance, chiral and electromagnetic current-conservation, etc. The Dyson-Schwinger 
equations (DSEs) provide such a framework. Experimental and theoretical studies are 
underway that will identify observable signals of M{p^) and thereby explain the most 
important mass-generating mechanism for visible matter in the Standard Model. 

There are many reasons why this is an exciting time in hadron physics. We have 
focused on one. Namely, through the DSEs, we are positioned to unify phenomena as 
apparently diverse as the: hadron spectrum; hadron elastic and transition form factors, 
from small- to large-Q^; and parton distribution functions. The key is an understanding 
of both the fundamental origin of nuclear mass and the far-reaching consequences of the 
mechanism responsible; namely, DCSB. These things might lead us to an explanation 
of confinement, the phenomenon that makes nonperturbative QCD the most interesting 
piece of the Standard Model. 
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